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ABSTRACT 

The evidence that stellar systems surrounding the Milky Way (MW) are distributed in 
a Vast Polar Structure (VPOS) may be observationally biased by satellites detected in 
surveys of the northern sky. The recent discoveries of more than a dozen new systems 
in the southern hemisphere thus constitute a critical test of the VPOS phenomenon. 
We report that the new objects are located close to the original VPOS, with half of 
the sample having offsets less than 20kpc. The positions of the new satellite galaxy 
candidates are so well aligned that the orientation of the revised best-fitting VPOS 
structure is preserved to within 9 degrees and the VPOS flattening is almost un¬ 
changed (31kpc height). Interestingly, the shortest distance of the VPOS plane from 
the MW center is now only 2.5 kpc, indicating that the new discoveries balance out 
the VPOS at the Galactic center. The vast majority of the MW satellites are thus 
consistent with sharing a similar orbital plane as the Magellanic Clouds, confirming a 
hypothesis proposed by Kunkel & Demers and Lynden-Bell almost 40 years ago. We 
predict the absolute proper motions of the new objects assuming they orbit within the 
VPOS. Independent of the VPOS results we also predict the velocity dispersions of 
the new systems under three distinct assumptions: that they (i) are dark-matter-free 
star clusters obeying Newtonian dynamics, (ii) are dwarf satellites lying on empirical 
scaling relations of galaxies in dark matter halos, and (iii) obey MOND. 

Key words: Galaxy: halo - galaxies: dwarf - galaxies: kinematics and dynamics - 
Local Group - Magellanic Clouds 


1 INTRODUCTION 

It has been known for almost 40 years now that the satel¬ 
lite galaxies of the Milky Way (MW) preferentially lie along 
a great circle in the sky which passes almost through the 
Galactic poles. The MW satellites are thus distributed in 
a flattened, polar plane. From the beginning it was noticed 
that this plane contains the Magellanic Clouds and follows 
the Magellanic stream. Hence, Kunkel V Demers (1976) and 
Lynden-Bell (1976) termed this phenomenon the “Magel¬ 
lanic plane”. As more MW satellites were discovered, e.g.. 
Sextans and those dwarf galaxies in the footprint of the 
Sloan Digital Sky Survey, they were found to lie close to 
this plane too (Kroupa, Theis V Boily 2005; Metz, Kroupa 
V Jerjen 2009). 

With the inclusion of other types of MW halo objects 
this plane was subsequently called the Vast Polar Struc¬ 
ture (VPOS) of the MW (Pawlowski, Pflamm-Altenburg 
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V Kroupa 2012). It was found that young halo globular 
clusters, which are hypothesised to have formed in external 
dwarf galaxies that were accreted into the Galactic potential 
well, follow the same polar distribution (Keller, Mackey V 
Da Costa 2012; Pawlowski, Pflamm-Altenburg V Kroupa 
2012). Even the stellar streams of some disrupting MW 
satellites (both galaxies and star clusters) have been found 
to align with the VPOS (Pawlowski, Pflamm-Altenburg V 
Kroupa 2012; Pawlowski V Kroupa 2014). The orbital direc¬ 
tions of the 11 classical satellites, deduced from their proper 
motions (PMs), further indicate that most satellites co-orbit 
in the plane. The VPOS is thus not only a spatial but also 
a rotating structure (Pawlowski V Kroupa 2013). 

Kroupa, Theis V Boily (2005) were the first to argue 
that the narrow, polar alignment of the 11 brightest “clas¬ 
sical” MW satellites is in conflict with the typical distri¬ 
bution of dark matter sub-halos in ACDM. Whether the 
positional alignment of these 11 satellite galaxies is prob¬ 
lematic for ACDM or not has been challenged since then 
(e.g. Deason et al. 2011; Libeskind et al. 2005, 2009; Wang, 
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Table 1. Compilation of observed properties of the newly detected MW satellite objects 


Object 

suggested names 

Type 

a 

s 

m - M 


Mv 

ri/2 

ri/2 

Ref. 




[°] 

[°] 

[mag] 

[kpc] 

[mag] 

[arcmin] 

[pc] 


Kim 1 

Kim 1 

GG 

332.92 

7.03 

16.5 ±0.1 

19.8 ±0.9 

0.3 ±0.5 

1.2 ±0.1 

6.9 ±0.6 

( 1 ) 

Ret II 

DES J0335.6-5403 

U 

53.92 

-54.05 

17.5 

32 

-3.6 ±0.1 

O o-|-1.0 
'^•^-0.6 

351^ 

(2) 


Reticulum 2 


53.93 

-54.05 

17.4 

30 

-2.7 ±0.1 

3.7 ±0.2 

32 ±1 

(3) 

Lae II 

Triangulum II 

U 

33.32 

36.18 

17.4 ±0.1 

30 

-1.8 ±0.5 

3 Q+i-i 

341^ 

(4) 

Tuc II 

DES J2251.2-5836 

UD 

343.06 

-58.57 

18.8 

58 

-3.9 ±0.2 

7.2 ±1.8 

120 ± 30 

(2) 


Tucana 2 


342.97 

-58.57 

19.2 

69 

-4.4 ±0.1 

9.9 ±1.4 

199 ± 28 

(3) 

Hor II 

Horologium H 

U 

49.13 

-54.14 

19.5 

78 

—2 


47± 10 

(5) 

Hor I 

DES J0255.4-5406 

u 

43.87 

-54.11 

19.7 

87 

-3.5 ±0.3 


60“^^® 

'^'^-30 

(2) 


Horologium 1 


43.88 

-54.12 

19.5 

79 

-3.4 ±0.1 

1.3 ±0.2 

30 ±3 

(3) 

Phe II 

DES J2339.9-5424 

u 

354.99 

-54.41 

19.9 

95 

-3.7 ±0.4 

1.20 ±0.6 


(2) 


Phoenix 2 


355.00 

-54.41 

19.6 

83 

-2.8 ±0.2 

1.1 ±0.2 

27 ±5 

(3) 

Eri III 

DES J0222.7-5217 

u 

35.69 

-52.28 

19.9 

95 

-2.4 ±0.6 

Q 4-t0.3 

lll« 

(2) 


Eridanus 3 


35.69 

-52.28 

19.7 

87 

-2.0 ±0.3 

0.7 ±0.3 

18 ±8 

(3) 

Kim 2 

Kim 2 

GG 

317.21 

-51.16 

20.1 ±0.1 

104.7 ±4.1 

-1.5 ±0.5 

0.42 ± 0.02 

12.8 ±0.6 

(6) 


DES J2108.8-5109 


317.20 

-51.16 

19.2 

69 

-2.2 ±0.5 

0.6 ±0.1 

12 ±2 

(2) 


Indus 1 


317.20 

-51.17 

20.0 

100 

-3.5 ±0.2 

1.4 ±0.4 

39± 11 

(3) 

Gru I 

Grus 1 

u 

344.18 

-50.16 

20.4 

120 

-3.4 ±0.3 

2.0 ±0.7 

70 ±23 

(3) 

Pic I 

DES J0443.8-5017 

u 

70.95 

-50.28 

20.5 

126 

-3.7 ±0.4 

-'■•^-0.6 

43 + 153 

(2) 


Pictoris 1 


70.95 

-50.28 

20.3 

114 

-3.1 ±0.3 

0.9 ±0.2 

31 ±7 

(3) 

Hyd II 

Hydra H 

UD 

185.44 

-31.98 

20.64 ±0.16 

134 ± 10 

-4.8 ±0.3 

1 7+0.3 
' -0.2 

68± 11 

(7) 

Peg III 

Pegasus HI 

UD 

336.10 

5.41 

21.56 ±0.20 

205 ± 20 

-4.1 ±0.5 

1.85 ±0.10 

110 ±6 

(8) 

Eri H 

DES J0344.3-4331 

UD 

56.09 

-43.53 

22.6 

330 

-7.4 ±0.1 

1.6 ±0.6 

155 ± 54 

(2) 


Eridanus 2 


56.09 

-43.53 

22.9 

380 

-6.6 ±0.1 

1.6 ±0.1 

172 ± 12 

(3) 


Properties of the newly discovered MW satellite objects collected from the literature, the reference from which the data of a given row 
has been collected is indicated in the last column (labelled “Ref.”). Type can be CC for confirmed star cluster, UD for unconfirmed dwarf 
galaxy candidate and U for unclassified (unconfirmed star cluster or ultra-faint dwarf galaxy candidate). The other parameters are a and 
5 for right ascension and declination of the object, m - M for its distance modulus, tq for the Heliocentric distance in kpc, Mv^ for the 
absolute v-band magnitude, and ri/2 f^e two-dimensional half light radius in both arc minutes and parsec. DES15 find Ret II and 
Erill to have a significant ellipticity e and report their half-light radii as measured along the semi-major axis. We report and use the 
circularized ri^2: which is obtained by multiplying the semi-major axis ellipticity with \/l — e. Uncertainties are given if they are specified 
in the original publication. 

Reference: (1): Kim & Jerjen (2015a), (2): DES15, (3): Koposov et al. (2015a), (4): Laevens et al. (2015), (5): Kim & Jerjen (2015b); (6): 
Kim et al. (2015b), (7): Martin et al. (2015), (8): Kim et al. (2015a). 


Frenk Sz Cooper 2013; Zentner et al. 2005). Taking into ac¬ 
count the limitation of detecting satellite galaxies close to 
the Galactic plane, the significance of the planar arrange¬ 
ment of the 11 satellites is between 99.4 and 99.9 per cent 
(Pawlowski et al. 2014; Pawlowski & McGaugh 2014a). Due 
to minor anisotropies present in dark matter sub-halo sys¬ 
tems around MW equivalents in simulations, it has been 
found that the VPOS flattening alone can be reproduced by 
0.5 to 6 per cent of ACDM systems (Wang, Frenk & Cooper 
2013; Pawlowski et al. 2014; Pawlowski & McGaugh 2014a). 
However, once the kinematic correlation is being taken into 
account, structures like the VPOS are extraordinarily rare 
in cosmological simulations (< 0.1 per cent, Pawlowski et al. 
2014; Pawlowski V McGaugh 2014a). Planar satellite ar¬ 
rangements sometimes occurring in the simulations tend to 
be transient features without aligned orbits (e.g. Gillet et al. 
2015). 

The significance of the VPOS can be tested with the 
help of additional satellite galaxies. In the data of the Sloan 
Digital Sky Survey (SDSS York et al. 2000) 15 faint and 


ultra-faint satellite galaxies have been discovered between 
2005 and 2010 (Willman et al. 2005a,b; Belokurov et al. 
2006, 2007, 2008, 2009, 2010; Sakamoto Sz Hasegawa 2006; 
Zucker et al. 2006b,a; Walsh, Jerjen V Willman 2007; Grill- 
mair 2009). These additional objects have supported the 
notion that the majority of MW satellites are part of a sin¬ 
gle halo structure: they independently define a polar plane 
which is closely aligned with the original disk of satellites 
(Kroupa et al. 2010). However, the VPOS results are based 
on the SDSS survey footprint, which covers the north Galac¬ 
tic pole region and thus the SDSS satellites may be expected 
to lie close to the VPOS. The SDSS survey nevertheless pro¬ 
vides additional constraints. Metz, Kroupa Sz Jerjen (2009) 
noted a deficit of SDSS-discovered satellite galaxies away 
from the satellite plane at large galactocentric distances. If 
the satellites discovered in the SDSS would be drawn from 
an isotropic distribution, the chance of them being oriented 
in a similarly narrow plane this closely aligned with that 
defined by the 11 classical satellites is small (Pawlowski in 
prep.). In addition, even though the SDSS survey area was 
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Figure 1. All-sky plot illustrating the positions of the MW satellites relative to the MW (inverted image in the background of the 
plot). Previously known satellite galaxies are plotted as filled black points, recently discovered satellite objects as white points and star 
clusters as filled black triangles. The velocities of the 11 classical satellites are indicated by the red arrows, which are the projected 3D 
velocity vectors determined from the measured line-of-sight velocities and proper motions as compiled in Pawlowski Sz Kroupa (2013). 
The dotted line indicates the orientation of the VPOS+new fit. The rms height of the VPOS is illustrated with two solid lines, which 
show the intersection of the VPOS+new of height 31 kpc with a sphere of 100 kpc radius. The majority of MW satellites fall close to the 
VPOS, as does the Magellanic Stream (dark patches). Note that the parallax effect due to the offset of the Sun from the Galactic center 
(and thus the best-fit VPOS plane) and the intrinsic thickness of the VPOS increase the scatter of the satellites around the VPOS great 
circle. The proper motions indicate that eight of the 11 classical satellites are consistent with co-orbiting in the VPOS, while Sculptor 
is counter-orbiting within the VPOS and Sagittarius is on an orbit perpendicular to both the MW and the VPOS, as also traced by 
the Sagittarius stream (contour lines indicate particle densities for the simulated stream from Law & Majewski 2010). The background 
image is by Nidever et al. (2010), NRAO/AUI/NSF and Meilinger, Leiden-Argentine-Bonn Survey, Parkes Observatory, Westerbork 
Observatory, Arecibo Observatory (see http://www.nrao.edu/pr/2010/magstream/). 


extended since data release 7, no new MW satellite galaxies 
were discovered outside of the VPOS (Pegasus III, recently 
discovered in SDSS DRIO data by Kim et al. 2015a, aligns 
well with the VPOS as will be shown later). 

Recently a number of new surveys began to operate, fa¬ 
cilitating a wider search for MW satellites beyond the SDSS 
footprint. The Pan-STARRSl survey covers three quarters 
of the entire sky {S > —30°), with much of this area be¬ 
ing far away from the VPOS and not covered by SDSS. 
Two new MW satellite objects have been discovered in this 
survey so far: the remote globular cluster PSO J174.0675- 
10.8774 (Laevens et al. 2014) also known as Crater (Be¬ 
lokurov et al. 2014), which happens to lie close to the VPOS 
(Pawlowski & Kroupa 2014), and the dwarf satellite Trian¬ 
gulum II (Laevens et al. 2015). 

In recent months, the Stromlo Milky Way Satellite Sur¬ 
vey (Jerjen 2010), the Dark Energy Survey (DES; The Dark 
Energy Survey Collaboration 2005) and the Survey of the 
Magellanic Stellar History (SMASH; Olsen et al. 2014; PI 


D. Nidever) have been revealing a series of new Milky Way 
companions: two ultra-faint star clusters Kim 1 (Kim & Jer¬ 
jen 2015a) and Kim 2 (Kim et al. 2015b), and two ultra-faint 
dwarf galaxies: Pegasus HI (Kim et al. 2015a) and Hydra H 
(Martin et al. 2015). Two studies of the DES-YlAl survey 
data announced the discovery of seven (The DES Collab¬ 
oration et al. 2015, , hereafter DES 15) respectively eight 
(Koposov et al. 2015a) new objects, not counting Kim 2 (or 
Indus I), which has been already found earlier (Kim et al. 
2015b). Eurthermore, Kim & Jerjen (2015b) have discov¬ 
ered an additional object. Horologium H, in the same data 
set which apparently was overlooked previously. One of the 
objects discovered in the DES, Erill, is most likely beyond 
the virial radius of the MW. Eundamental parameters for 
these new objects are compiled in Table 1. Eig. 1 shows 
their positions in relation to the MW, the previously known 
satellites and the Magellanic Stream in an all-sky plot. 

Many of the new objects are candidate satellite galax¬ 
ies although the distinction between star clusters and ultra- 
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faint satellite galaxies based on their size and luminosity 
alone is becoming increasingly blurry. We therefore refrain 
from choosing a preferred name, but instead list all sug¬ 
gested names in Table 1. We will identify the objects by the 
abbreviations given in the hrst column, consisting of three 
letters and a Roman numeral, except Kim 1 & 2, which have 
been unambiguously classified as star clusters. 

The alignment of the satellite galaxies in a common 
structure indicate that they might be dynamically asso¬ 
ciated, sharing similar orbits. This assumption allows us 
to predict their PMs (Lynden-Bell & Lynden-Bell 1995; 
Pawlowski & Kroupa 2013), measurement of which then pro¬ 
vides a crucial test of whether the assumed association is 
real. The existence of the VPOS and the preferential align¬ 
ment of streams with this structure indicates such a dynam¬ 
ical association of its constituents. The fact that the 11 clas¬ 
sical MW satellites already have PM measurements which 
show that most of them indeed (co-)orbit in the VPOS fur¬ 
ther supports this interpretation. The PMs of the new MW 
satellite objects can thus be predicted in an entirely em¬ 
pirical way (Pawlowski V Kroupa 2013). This prediction is 
based solely on the current spatial distribution of the MW 
satellites and does not require an assumption of a specific 
MW potential or underlying type of dynamics. Confidence in 
the method can be drawn from the finding that more precise 
PM measurements tend to agree better with the predicted 
PMs (Pawlowski V Kroupa 2013). 

The observed distribution of stellar light (and thus 
stellar mass) in the newly discovered objects can be used 
to predict the velocity dispersions of the newly discovered 
objects. The results depend strongly on which dynamical 
model (Newtonian dynamics or Modified Newtonian Dy¬ 
namics, MOND, Milgrom 1983) or which dark matter halo 
scaling relation is assumed (McGaugh et al. 2007; Walker 
et al. 2009). We will provide predictions for all these differ¬ 
ent cases. 

This paper is structured as follows. In Sect. 2 we will 
compare the positions of the new satellite objects with the 
VPOS defined by the previously-known MW satellites be¬ 
fore measuring the effect of adding the new objects to the 
VPOS plane ht. Our predictions for the proper motions of 
the objects are presented in Sect. 3, and the predicted veloc¬ 
ity dispersions in Sect. 4. We end with concluding remarks 
in Sect. 5. 


2 THE VPOS 

If the VPOS is a real structure, then it makes a prediction 
for the probable locations of new satellites. Here we compare 
the positions of the new satellite objects with the VPOS as 
dehned by the previously-known MW satellite galaxies (sect. 
2.1). We then determine how the VPOS orientation changes 
when the new objects are included in the plane fit. 

For Kiml, Lae II, Horll, HydII, Trill and Peg III we 
will use the positions and distances reported in their re¬ 
spective discovery papers, as listed in Table 1. For those 
satellite objects discovered by both DES15 and Koposov 
et al. (2015a) we will use the average of the two distances 
estimates. The one exception is the star cluster Kim 2, for 
which more reliable observational data was obtained in the 
discovery paper by Kim et al. 2015b, such that we use their 


distance measurement. We found that our results are statis¬ 
tically robust against either using the DES15 or the Koposov 
et al. (2015a) distances. 

2.1 Comparison to the known VPOS 

As a first step, we have compared the positions of the new 
MW satellite objects with the VPOS plane orientation given 
by (i) the 27 previously known MW satellite galaxies (here¬ 
after VPOSall) and (ii) the ht to 24 of those galaxies, ex¬ 
cluding the three outliers Hercules, Ursa Major and Leo I 
(hereafter VPOS-3). The respective plane parameters from 
Pawlowski, Kroupa V Jerjen 2013 are reproduced in the sec¬ 
ond and third column of Table 2. As already mentioned by 
Kim V Jerjen (2015a), Kim et al. (2015b), DES15, Martin 
et al. (2015) and Kim et al. (2015a), the new objects lie close 
to the VPOS. Table 3 lists their offsets from the VPOSall 
and VPOS-3. The average offset of the 13 new objects from 
the VPOSall plane is only 31 kpc, the median is 21 kpc. 

Eight objects are particularly close to the best-ht 
planes. The seven satellite objects Kiml, Retll, Horl, 
Horll, Erilll, Pic I and Peg HI have offsets of < 20 kpc from 
the VPOSall, and all but Eri H have offsets of only 11 kpc 
or less from the VPOS-3. Considering Eri IPs large distance 
of 350 ± 30 kpc from the MW, it is worth noticing that the 
two distance ratios dyposaii/dMW and dvpos-s/dMW of less 
than 14 percent are remarkably small. 

The largest offsets are found for the candidate objects 
Gru I, Hyd H and the star cluster Kim 2. Koposov et al. 
(2015a) note that Gru I is located close to a CCD chip gap in 
the DES survey, such that its parameters are more uncertain 
than those of the other candidates. A smaller distance from 
the MW would reduce the offset. Eurthermore, in particular 
the two dwarf galaxy candidates Gru I and Hyd H both 
have large Galactocentric distances, such that despite their 
relatively large perpendicular offset from the planes their 
position vectors are only about 30° inclined relative to the 
planes. In summary, most of the newly reported satellite 
galaxies and star clusters are qualitatively consistent with 
the original VPOS. 

2.2 Influence of the DES footprint 

The region covered by the DES-YlAl search area is close 
to the Magellanic Clouds, and thus not far from the VPOS 
plane. It might therefore not be surprising that the satellite 
objects discovered by this survey are close to the VPOS, 
as argued by DES 15. To quantify how much of an align¬ 
ment is expected due to the survey footprint, and thus to 
test whether it causes the alignment, we have created 10000 
realizations of isotropic DES object positions. To construct 
these, the angular positions of the ten DES objects have been 
randomly selected from an isotropic distribution around the 
Galactic center, while their Galactocentric distances have 
been preserved. Only positions lying within the region cov¬ 
ered by the already observed part of the DES footprint were 
accepted. Positions outside of this region were randomized 
again until they were within the survey footprint. This guar¬ 
antees that each randomized realisation contains ten objects 
which share the same radial distribution as the observed ob¬ 
jects. We have then determined the offsets from the VPOS¬ 
all plane for each position, and determined the mean and 
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Table 2. New plane fit parameters (all objects except confirmed star clusters) 


Name VPOSall VPOS-3 VPOS+new VPOS+new-4 VPOSsouth VPOSnorth 



^MW [kpc] 

A [kpc] 

c/a 

b/a 

A^members 

Outliers excluded? 
Includes new objects? 


155.6\ /169.5\ 

-3.3y V-2.8y 


7.9 

10.4 

29.3 

19.9 

0.301 

0.209 

0.576 

0.536 

27 

24 

no 

yes 

no 

no 


d64.0\ 

/169.4' 

-6.9 ) 

V-6-1. 

2.5 

5.4 

30.9 

21.3 

0.313 

0.224 

0.579 

0.566 

38 

34 

no 

yes 

yes 

yes 


169.3\ /157.7\ 

-7.2y V-8.9y 


3.6 

3.0 

22.3 

36.4 

0.356 

0.475 

0.619 

0.623 

19 

19 

no 

no 

yes 

yes 


Parameters of the plane fits: 

n: The direction of the normal vector (minor axis) of the best-fit plane in Galactic longitude I and latitude b. 

Dyiw'- offset of the planes from the center of the MW. 

A: RMS height from the best-fit plane of the MW satellite objects included in the fit. 

c/a and b/a: short- and intermediate-to-long axis ratios, determined from the RMS heights in the directions of the three axes. 
A^members- Number of objects used for the fits. 


median offsets for each set of ten randomized DES satellite 
objects. 

Figure 2 shows the resulting distribution of the mean 
and the median offsets of these 10000 realisations. The 
ten observed objects have an average offset of 28 kpc and 
a median offset of 18 kpc. If they were drawn from an 
isotropic distribution, the expected (average) mean and 
median would be considerably (about 50 per cent) larger: 
41 kpc and 31 kpc, respectively. Of all randomized realisa¬ 
tions, 91 per cent have a larger mean offset, and 96 per cent 
have a larger median offset than the observed values. This in¬ 
dicates that the observed alignment is indeed stronger than 
expected from the current survey footprint alone, but the 
significance of this conclusion is not extremely high. This 
might change once the full DES footprint is covered, in par¬ 
ticular if at larger distances from the VPOS fewer satellites 
are discovered. 


2.3 Effects on the best-fit VPOS 

We now update the plane fits by including the new objects. 
Since we use the same fitting routine and adopt the same 
parameters as in Pawlowski, Kroupa V Jerjen (2013), we 
refer the reader to that paper for further information on the 
method. 

To be consistent with the previous plane fitting analysis 
that focussed on satellite galaxies, in the following we add 
only those 11 new objects to the plane fit sample which 
are likely but still unconfirmed satellite galaxies. The two 
objects Kim 1 V 2, which have already been identified as star 
clusters through follow-up observations are not included. We 
also exclude Erill from the fit because its large distance 
estimate places it outside of the virial radius of the MW 
(?^ 250kpc), thus should not be considered a satellite galaxy 
of the Milky Way. 

For the full VPOS sample (called VPOS+new in Table 
2) we use all 27 previously identified satellites plus the 11 
new satellite objects. Our analysis shows that even though 
the sample size used for the fit has increased by 40 percent, 
and the number of objects in the southern Galactic hemi¬ 
sphere doubled, the VPOS parameters remain essentially 


unchanged. The direction of the plane normal n changes by 
only about 9°, to larger Galactic longitude and slightly lower 
Galactic latitude. This normal is closely aligned (within 9°) 
with the normal direction of the ‘classical’ disk of satel¬ 
lites dehned by the 11 brightest MW satellites, which points 
to (/,6) = (157.3°,—12.7°). Because these are the bright¬ 
est MW satellites, their distribution should be the least 
affected by biases due to uneven sky coverages, but low- 
number statistics are a concern. 

Our plane ht routine does not require the htted plane 
to pass through the MW center. This extra freedom acts 
as a consistency check for whether a found plane can be 
dynamically stable^. It is therefore arguably the most in¬ 
teresting hnding that the new offset from the MW center 
Dmw, is reduced from 7.9 to only 2.5 kpc. This might indi¬ 
cate that it needed more satellites in the southern Galactic 
hemisphere to ‘balance’ VPOS out in the MW center. De¬ 
spite the overall shift to smaller MW offsets, both measures 
of the plane thickness, the root-mean-square height A and 
the short-to-long axis ratio c/a, remain almost unchanged 
(they both increase by a few per cent) if the additional satel¬ 
lite galaxy candidates are included in the ht. The resulting 
values A = 30.9 kpc and c/a = 0.313 are both at least 50 per 
cent smaller than those expected for the same number of 
satellites objects drawn from isotropic distributions or from 
sub-halo distributions in cosmological simulations of MW 
equivalents (see e.g. hg. 3 of Pawlowski V McGaugh 2014a). 

The distances of the satellites from this new best-ht 
plane are also compiled in Table 3. Not surprisingly, the new 
objects have smaller offsets from the VPOS plane if they are 
included in the plane ht. Since the plane orientation has not 
changed dramatically compared to the VPOSall, the same 
objects tend to be closest to the new ht, too. 

In addition to htting a plane to ah 27 known plus the 11 
new satellite candidates, we have also constructed a sample 
analogous to the VPOS-3 by excluding the four outliers with 
offsets of more than 50 kpc from the VPOS+New ht. These 

^ If a satellite plane has a large offset from its host’s center (such 
as Plane 2 in Shaya & Tally 2013), the satellites can not orbit 
within the plane and the arrangement must be a transient feature. 
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Figure 2. Histogram (upper panel) and cumulative distribution 
(lower panel) of mean and median offsets for 10000 randomized 
distributions of the 10 stellar systems within the DES-YlAl foot¬ 
print, as discussed in Sect. 2.2. Even though the DES footprint 
is close to the VPOS, the observed mean and median (thick solid 
blue and thinner dashed red histograms, respectively) offsets of 
the 10 objects (vertical dashed and dotted lines, respectively) are 
smaller than expected for isotropically distributed objects. 

are the three previously-identified outliers Hercules, Ursa 
Major I, and Leo I, as well as the newly discovered object 
Hyd II. We note that Leo I has a proper motion indicating 
that it does not orbit within the VPOS. It might not even be 
a true MW satellite given its large Galactocentric distance 
and high velocity. This VPOS+New-4 has a very similar 
orientation as the VPOS+New (inclined by about 6°) but 
is significantly thinner (A = 21.3 kpc, cja — 0.224). It is 
almost identical in orientation and thickness to the VPOS- 
3, but also has a smaller offset from the MW center. 

As reported before (Pawlowski, Kroupa h Jerjen 2013), 


VPOS+new (edge-on) 



VPOS-i-new (face-on) 



[kpc] 

Figure 3. Distribution of satellite objects in Cartesian coordi¬ 
nates around the MW. The upper panel shows the VPOS+new 
fit edge-on (black solid line, the dashed lines indicate the rms 
height), the lower panel shows a view rotated by 90°, in which 
the VPOS is oriented approximately face-on. In this view, the 
MW satellites with measured proper motions preferentially or¬ 
bit in the clockwise direction. The 11 brightest (classical) MW 
satellite galaxies are plotted as yellow dots, the fainter satellite 
galaxies as smaller green dots and globular clusters classified as 
young halo objects as blue squares. New objects confirmed to be 
star clusters (PSO J174.0675-10.8774/Crater in the north, Kim 1 
& 2 in the south) are plotted as lighter blue diamonds, all other 
new objects as bright-green stars. The red lines in the center in¬ 
dicate the position and orientation of streams in the MW halo. 
They preferentially align with the VPOS, but are mostly con¬ 
fined to the innermost regions of the satellite distribution. Both 
plots are centred on the MW (cyan line) which is seen edge-on. 
The grey wedges indicate the region (±12° around the MW disk, 
where satellite galaxies might be obscured by the Galaxy. 
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the VPOS plane is oriented almost exactly like the cur¬ 
rent orbital plane of the LMC, which has an orbital pole 
(direction of angular momentum) pointing towards {l,b) = 
(175.4°,—5.7°) (Pawlowski V Kroupa 2013), less than 12° 
inclined with respect to the VPOS+new and only 6° with 
respect to the VPOS+New-4. This supports the almost 40 
years old notion by Kunkel V Demers (1976) and Lynden- 
Bell (1976) that most MW satellites seem to share the same 
orbital plane as the Magellanic Clouds. 

Fig. 3 shows the distribution of known satellite ob¬ 
jects around the MW from both an edge-on and face-on 
views of the updated VPOS+new. Newly discovered satel¬ 
lite objects are plotted as light blue diamonds if they have 
been conhrmed to be star clusters, or as bright green stars. 
It is obvious that most of the recently discovered satellite 
galaxy candidates lie within the VPOS plane at intermedi¬ 
ate distances from the MW. The majority populate a region 
where previously-known classical MW satellites (yellow dots 
in the figure) are situated^. This is illustrated by the face-on 
VPOS+New view in the lower panel of Fig. 3, in which the 
majority of new discoveries populate the lower (southern) 
left quadrant of the VPOS plane. 

This panel also illustrates how isolated the new MW 
satellite Pegasus III is. As pointed out by Kim et al. (2015a), 
its only known close neighbour at a distance of about 37 kpc 
is Pisces II. The next-closest known MW satellites are found 
at distances of 170 kpc or more. The whole lower (southern) 
right quadrant of the VPOS in the lower panel of Fig. 3 
appears to be surprisingly devoid of known MW satellites. 
None of the 11 classical MW satellites falls into this region 
at Galactocentric distances beyond about 30 kpc (the clos¬ 
est one, Sagittarius, lies at the edge of this quadrant but is 
known to orbit perpendicular to the VPOS). Interestingly, 
part of this quadrant was surveyed by the SEGUE extension 
of the SDSS survey, within which Pisces II, the only previ¬ 
ously known faint MW satellite galaxies in this region, was 
found. The next-closest known structure, also found in the 
data of the SDSS survey, is the Pisces overdensity (Watkins 
et al. 2009) at a distance of 130 kpc from Pegasus III (green 
triangle in Eig. 3). The overdensity might be the remnant of 
a dwarf galaxy or star cluster currently being disrupted, and 
is also aligned with the Magellanic Stream and the VPOS. 


2.4 Dividing the VPOS: north-south differences? 

Eor a long time, the majority of known MW satellites galax¬ 
ies lay in the Galactic north because the SDSS initially con¬ 
centrated on that part of the sky. With the discovery of 
new satellite galaxy candidates, the numbers of southern ob¬ 
jects approaches that of the northern ones (19 each, if all of 
the new candidate objects turn out to be bona-hde satellite 
galaxies). The comparable sample size allows us to fit the 
VPOS for the southern and northern satellite populations 
separately. Such a split might reveal whether the VPOS is 
systematically tilted or bent, for example due to precession 
of the satellites on their orbits or because of non-symmetric 
infall of the VPOS-satellites. Due to the Galactic disc which 
possibly obscures satellites along the Galactic equator, split- 


^ The LMC and SMC, Carina, Fornax and Sculptor 


ting the satellites into a northern and southern sample is the 
most sensible separation. 

The resulting plane-fit parameters for the two sub¬ 
samples are similar to each other and to the VPOS+New 
plane fitted to all objects (see Table 2). Both the northern 
and the southern VPOS plane fits are inclined by only 6° 
relative to the total VPOS+New plane fit. They both are 
almost polar (the normal vectors have Galactic latitudes of 
b = —7° and —9°, for the southern and the northern satel¬ 
lites, respectively), have a very similar orientation (inclined 
by less than 12°), have essentially the same small offset 
from the Galactic center (Dmw ~ 3 kpc) and are similarly 
thin (Arms = 22 versus 36kpc). However, while the thinner, 
southern plane shares its orientation with the orbital plane 
of the LMG (and the VPOS-3 plane fits) to within 6°, the 
northern one is more inclined (18°). 

Interestingly, the orientation and thickness of the north¬ 
ern plane is strongly affected by the four outliers, which are 
all northern objects. If they are excluded from the north¬ 
ern VPOS ht, the resulting plane parameters are extremely 
similar to those of the southern fit: the normal to the best- 
fit plane points to (/,5) = (168.5°,—3.3°), the root-mean- 
square height of the plane is Arms = 20 kpc, but the offset 
from the MW center is Dmw = 11 kpc. Even without exclud¬ 
ing the outliers from the northern satellite sample, the plane 
parameters of the northern (VPOSnorth) and the southern 
(VPOSsouth) fits agree well. We therefore conclude that the 
current data suggest that the northern and the southern part 
of the VPOS have similar properties and orientations. This 
might also indicate that observational biases - which are dif¬ 
ferent for the northern and southern hemispheres - do not 
affect the determined plane orientations significantly. 

2.5 The distant object Eri II 

The two current distance estimates for Eri II are 330 
(DES15) or 380kpc (Koposov et al. 2015a). These place the 
object barely outside of the virial radius expected for the 
MW (about 250 to 300kpc). If its position just beyond the 
MW halo is confirmed, this would make Eri II the closest 
known non-satellite dwarf galaxy in the Local Group, closer 
than the neighbours Phoenix and Leo T at about 420 kpc 
distance. 

Eri II aligns well with the VPOS plane, even though it 
was excluded from the fit due to its large distance. In this 
regard it might be interesting to check where Eri II lies with 
respect to the dominating plane of non-satellite galaxies in 
the LG, termed LG plane 1 (LGPl, Pawlowski, Kroupa Sz 
Jerjen 2013; Pawlowski V McGaugh 2014b). LGPl consists 
of about 10 LG dwarf galaxies which have distances of more 
than 300 kpc from both the MW and M31. It is an extremely 
narrow plane, with root-mean-square height of about 50 kpc 
and a maximum diameter of 2Mpc, which stretches from 
M31 to the MW and beyond. Approximately in Eri IPs di¬ 
rection LGPl has its closest approach to the MW, of about 
180 kpc. Depending on its exact position, Eri II has an offset 
of about 100 kpc from LGPl. It therefore may potentially be 
a member of this structure. Assuming Eri II to be part of 
and perfectly aligned with LGPl, we would have predicted 
it to be closer than the current estimates by about 50 to 
100 kpc. Measuring Eri IPs line-of-sight velocity might shed 
more light on whether an association with LGPl is likely. 
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Table 3. Galactocentric distances, Cartesian positions and offsets of MW satellites from the VPOS plane fits (in kpc) 


Name 

dmw 

X 

y 

z 

C^VPOSall 

C^VPOS-3 

C^VPOS+New 

C^VPOS+New-4 

^VPOSsouth 

C^VPOSnorth 

The Galaxy 

0.0 

0.0 

0.0 

0.0 

8 

10 

3 

5 

4 

3 

Ganis Major 

13.4 

-II.9 

-6.2 

-I.O 

1 

1 

7 

5 

7 

II 

Sagittarius dSph 

18.4 

I7.I 

2.5 

-6.4 

23 

28 

19 

22 

21 

12 

Segue (I) 

27.9 

-19.4 

-9.5 

17.7 

5 

7 

12 

11 

12 

15 

Ursa Major II 

38.0 

-30.6 

II.6 

19.2 

26 

23 

30 

27 

29 

35 

Bootes II 

39.5 

6.6 

-1.7 

38.9 

17 

19 

14 

17 

14 

10 

Segue II 

40.8 

-31.8 

13.9 

-21.4 

29 

26 

35 

32 

33 

41 

Willman I 

42.9 

-27.7 

7.6 

31.8 

22 

20 

25 

22 

25 

29 

Goma Berenices 

44.9 

-10.6 

-4.3 

43.4 

3 

3 

1 

1 

2 

4 

Bootes III 

45.8 

1.3 

6.9 

45.3 

9 

13 

7 

11 

8 

2 

LMG 

50.0 

-0.6 

-41.8 

-27.5 

24 

16 

11 

10 

9 

10 

SMG 

61.2 

16.5 

-38.5 

-44.7 

38 

33 

25 

25 

24 

23 

Bootes (I) 

64.0 

14.8 

-0.8 

62.2 

26 

29 

25 

27 

25 

20 

Draco 

75.9 

-4.4 

62.2 

43.2 

21 

4 

15 

6 

9 

27 

Ursa Minor 

77.8 

-22.2 

52.0 

53.5 

33 

19 

28 

21 

24 

38 

Sculptor 

86.0 

-5.2 

-9.8 

-85.3 

2 

3 

10 

7 

8 

16 

Sextans (I) 

89.0 

-36.7 

-56.9 

57.8 

2 

13 

9 

14 

16 

5 

Ursa Major (I) 

I0I.6 

-61.1 

19.8 

78.7 

53 

51 

54 

51 

54 

57 

Garina 

106.8 

-25.1 

-95.9 

-39.8 

25 

1 

2 

6 

6 

9 

Hercules 

126.1 

84.1 

50.7 

79.1 

71 

94 

83 

93 

89 

68 

Fornax 

149.3 

-41.3 

-51.0 

-I34.I 

17 

30 

40 

42 

41 

40 

Leo IV 

154.8 

-I5.I 

-84.8 

128.6 

39 

18 

29 

21 

18 

37 

Ganes Venatici H 

160.6 

-16.5 

18.6 

158.7 

6 

2 

1 

3 

5 

4 

Leo V 

178.6 

-21.5 

-91.9 

I5I.7 

37 

14 

28 

18 

16 

38 

Pisces H 

I8I.I 

14.9 

I2I.7 

-133.3 

38 

4 

34 

17 

19 

57 

Ganes Venatici (I) 

217.5 

2.1 

37.0 

214.3 

6 

17 

20 

24 

19 

14 

Leo H 

235.9 

-77.3 

-58.3 

215.2 

26 

46 

30 

38 

42 

20 

Leo I 

257.4 

-123.6 

-II9.3 

I9I.7 

45 

83 

60 

76 

79 

41 

Kim I 

I9.I 

-2.7 

14.4 

-12.3 

2 

4 

6 

I 

3 

13 

Ret H 

33.0 

-9.7 

-20.4 

-24.1 

6 

3 

4 

3 

4 

7 

Lae 2 (Tri H) 

36.6 

-29.8 

17.4 

-12.2 

27 

23 

32 

28 

30 

39 

Tuc H 

59.2 

24.4 

-20.4 

-49.9 

36 

35 

26 

28 

27 

20 

Hor H 

80.0 

-9.6 

-48.7 

-62.7 

16 

7 

1 

1 

2 

4 

Hor I 

83.5 

-7.2 

-48.0 

-67.9 

18 

9 

2 

0 

1 

3 

Phe H 

88.1 

28.7 

-27.2 

-78.8 

41 

40 

28 

30 

29 

23 

Eri HI 

91.2 

-4.3 

-46.0 

-78.7 

19 

II 

3 

1 

1 

4 

Kim H (Ind I) 

98.9 

67.5 

-17.3 

-70.2 

72 

76 

63 

67 

66 

55 

Grus I 

II6.4 

50.6 

-23.0 

-102.2 

57 

59 

45 

48 

48 

38 

Pic I 

I2I.9 

-28.1 

-88.2 

-79.2 

15 

5 

8 

14 

14 

6 

Hydra H 

129.0 

40.8 

-102.4 

66.8 

92 

72 

79 

71 

70 

85 

Pegasus HI 

203.1 

44.3 

143.5 

-136.8 

21 

21 

13 

8 

6 

40 

Eri H 

365.0 

-86.2 

-2II.4 

-284.7 

4 

50 

51 

69 

66 

35 


MW satellite distances from the MW (dMw)^ Galactocentric Cartesian x, y, and z positions and offsets from the different best-fit VPOS 
planes: all 27 satellite galaxies considered in Pawlowski, Kroupa & Jerjen (2013) (VPOSall), the same sample but excluding three outliers 
(VPOS-3) and the fit to the 27 known and confirmed MW satellite galaxies plus all new objects except those which have been identified as 
star clusters (VPOS+new, i.e. the VPOSall sample plus Ret II, Lae 2, TucII, Horl and II, Phell, Erilll, Grul, Pic I, HydII and Peg III). 
For those of the new objects discovered by both B15 and K15 we assume that their Heliocentric distances are the average of the two 
distance estimates when fitting the plane. If an object was included in the plane fit its respective offset if printed in boldface. 


All LGPl dwarf galaxies have line-of-sight velocities relative 
to the local standard of rest which follow the Magellanic 
Stream velocity at their positions (see Fig. 6). If Eri II fol¬ 
lows the same trend, its velocity should be in the range of 
50 to 250 kms“^. 

The spatial position of Eri II thus places it at the inter¬ 
section of LGPI and the VPOS, a very interesting position 
holding the promise that Eri II might provide important 
clues on if and how the VPOS is connected to the larger- 
scale structure of dwarf galaxies in the LG. 


3 PREDICTED PROPER MOTIONS 

Out of the II MW satellite galaxies for which proper mo¬ 
tions (PMs) have been measured, eight are consistent with 
co-orbiting in the same sense within the VPOS, while Sculp¬ 
tor also orbits within the VPOS but in the opposite di¬ 
rection (Metz, Kroupa Libeskind 2008; Pawlowski 
Kroupa 2013). This orbital alignment finds further sup¬ 
port in the alignment of several streams in the MW halo 
with the VPOS (Pawlowski, Pflamm-Altenburg Kroupa 
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Figure 4. Predicted proper motions of the new MW satellite objects assuming that they move within the VPOS-3 (which is closely 
aligned with the orbital plane of the LMC). The map illustrates the cosine of the angle 6 between the VPOS-3 plane (the plane fitted to 
all confirmed satellite galaxies except three outliers, see Pawlowski, Kroupa & Jerjen 2013 for more details) and the current orbital plane 
which would result from the PM for each combination of the components [i(x cos (5 and 115 . The radial grey contour lines illustrate cos^ 
in steps of 0.2. For cos 6 > 0.8 an object’s orbital plane is inclined by less than 37° with the VPOS-3, and can therefore be considered 
to be co-orbiting (best co-orbiting alignment is marked with a thick magenta line). For cos 6 < —0.8 it would be counter-orbiting (best 
alignment marked with a thin magenta line). The green circular contours indicate the lower limits on the absolute speed of the objects 
relative to the MW in kms“^. Since their kinematics have not yet been measured, this assumes that they have zero line-of-sigh velocity 
relative to the Galactic standard of rest, resulting in the Heliocentric line-of-sight velocities as listed in Table 4. 


© 2015 RAS, MNRAS 000 , 1-17 











































































10 Pawlowski, McGaugh & Jerjen 


Table 4. Predicted PMs of MW satellite objects assuming they orbit within the VPOS-3 


Name 


[kpc] 


min(uios) 

[kms“^] 


0 : 


predicted 

VPOS-3 


'’min 5 "dmax] 

[kms“^] 


lloc COS 5 
[mas yr“ 


fia. cos 5 
[mas yr“ 


Kim 1 

68.5 

-38.4 

20 

Ret II 

266.3 

-49.7 

32 

Tuc II 

328.1 

-52.3 

63 

Hor II 

268.5 

-52.2 

79 

Hor I 

271.4 

-54.7 

83 

Phe II 

323.7 

-59.7 

89 

Eri III 

275.0 

-59.6 

91 

Kim 2 

347.2 

-42.1 

105 

Gru I 

338.7 

-58.2 

120 

Hyd II 

295.6 

30.5 

132 

Pic I 

257.4 

-41.2 

120 

Peg III 

69.8 

-41.8 

205 

Eri II 

249.8 

-51.6 

363 


-182 

17.8 

[50, 595] 

168 

11.9 

[50, 518] 

81 

25.3 

[50, 455] 

160 

2.4 

[50, 429] 

151 

1.0 

[50, 421] 

77 

19.5 

[50, 423] 

132 

0.3 

[50, 409] 

38 

41.9 

[50, 390] 

49 

24.9 

[50, 369] 

187 

28.9 

[50, 361] 

191 

7.1 

[50, 361] 

-174 

2.9 

[50, 273] 

154 

9.3 

[50, 105] 


[+0.92, 

[-2.17, 

/[+1.43, 

V[-0.66, 

/[+0.33, 

Vi-0.90, 

/[+0.57, 

U-O-31, 

[+0.54, 

[-0.33, 

/[+0.31, 

l^[-0.59, 

[+0.49, 

[-0.37, 

/[+0.03, 

l^[-0.60, 

[+0.19, 

[-0.49, 

[-0.18, 

[-0.09, 

[+0.35, 

[- 0 . 11 , 

[+ 0 . 10 , 

[- 0 . 22 , 

[+ 0 . 12 , 

[-0.07, 


+ 4.55]\ 

- 6.72]J 
+ 4.41]\ 

- 1-61] j 

+ 0.74]\ 

- 2.20]J 

+ 1.50]\ 

- 0.68]y/ 

+ 1.39]\ 

- 0.74]^ 
+ 0.73]\ 

- 1.35]yl 
+ 1.19]\ 
-0.81]^ 
+ 0.08]\ 

- l-28]y' 
+ 0.41]\ 

- l.OOjyl 
+ 0.03]\ 
+ 0.36]^ 
+ 0.90]\ 
-o.iejyl 
+ 0.24]\ 
-0.40]) 
+ 0.15]\ 

-0.08]y) 


[+0.26, 

[-1.34, 

[+0.80, 

[-0.46, 

[+0.23, 

[-0.59, 

[+0.32, 

[- 0 . 21 , 

[+0.32, 

[- 0 . 22 , 

[+ 0 . 20 , 

[-0.39, 

[+0.29, 

[-0.24, 

[+ 0 . 01 , 

[-0.40, 

[+ 0 . 12 , 

[-0.32, 

[-0.25, 

[-0.24, 

[+0.18, 

[- 0 . 10 , 

[+0.04, 

[-0.13, 

[+0.07, 

-0.06, 


- 3.36]\ 
+ 3.21]J 

- 2.18]\ 
+ 0.48]^ 
- 0.18]\ 
+ 0.71]^ 
- 0.61]\ 
+ 0.17]^ 

- 0.53]\ 
+ 0.19]j 

- 0 . 21 ]\ 

+ 0.37]^ 

- 0.41]\ 

+ 0 . 20 ]^ 

- 0.04]\ 
+ 0.29]j 
- 0 . 10 ]\ 
+ 0.19]^ 

- 0.46]\ 

- 0.69] j 

- 0.36]\ 

- 0.05] j 
- 0 . 10 ]\ 
+ 0.05]^ 
+ 0.04]\ 

- 0.05]y) 


I, b, tq: Heliocentric position of the MW satellite object in Galactic longitude and latitude and radius. 

min(r;ios): heliocentric line-of-sight velocity that minimizes the line-of-sight component of its Galactocentric velocity (i.e. the negative of 
the Galactocentric solar velocity component in the direction of the object). 

^VPO+4 smallest-possible inclination between the orbital plane of the satellite galaxy and the VPOS-3 plane, defined by the angle 
between the Galactocentric position and the plane normal. 

Umax: the range of PMs is constrained by adopting this range of minimum and maximum absolute speeds for the objects. The 
minimum is set to 50kms“^, the maximum determined by requiring the object to be approximately bound to the MW (see Pawlowski & 
Kroupa 2013 for details). 


fia cos 5 
IJ'6 


and 


fia cos S 
IJ-5 


predicted PM range if co- or counter-orbiting. 


2012) , most prominently the Magellanic Stream. That satel¬ 
lite galaxy planes appear to be rotating structures is also 
supported by the M31 satellite plane which shows coher¬ 
ent line-of-sight velocities indicative of rotation (Ibata et al. 

2013) . In addition, the unexpected anti-correlation of veloc¬ 
ities for diametrically opposite satellite galaxy pairs in the 
SDSS indicates that co-orbiting satellite galaxy planes might 
be ubiquitous (Ibata et al. 2014). 

Since most of the newly discovered objects are closely 
aligned with the VPOS-3 plane and thus with the average 
orbital plane of the co-orbiting MW satellite galaxies, we 
can empirically predict the PM of the new MW satellite ob¬ 
jects by using the method presented in Pawlowski & Kroupa 
(2013). This method assumes that the object it either co- or 
counter-orbiting in the plane defined by the current satel¬ 
lite positions. Because no line-of-sight velocity has yet been 
measured for most of the new objects, we assume that they 
have minimum (i.e. zero) Galactocentric line-of-sight veloc¬ 
ities min(uios), calculated as the negative of the component 
of the Galactocentric Solar motion in the direction of the 


object. The exact line-of-sight velocity of an object does not 
affect the orientation of the predicted PM, but only con¬ 
strains the maximum and minimum predicted PM to val¬ 
ues such that it remains bound to the MW. By assuming a 
minimum Galactocentric line-of-sight velocity the predicted 
PM range will be maximal. Once line-of-sight velocities have 
been measured they can be used to put tighter constraints 
on the expected range in PMs. 

To make the prediction compatible with those presented 
in Pawlowski & Kroupa (2013) and Pawlowski & Kroupa 
(2014), we adopt the same values for the circular velocity of 
the MW, the solar motion with respect to the local standard 
of rest and the Galactocentric distance of the Sun. This is 
also one reason why the orientation of the VPOS-3 was used 
in this prediction instead of the new fit. Another reason is 
that it is not yet known with certainty which of the new ob¬ 
jects are satellite galaxies and which are star clusters that 
should not be part of the VPOS fit. However, since the mea¬ 
sured VPOS orientation changes only minimally if the new 
objects are included in the plane fits, this choice does not 
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strongly affect the predicted PMs anyway. We again use the 
average of the distances reported by K15 and B15 for those 
objects discovered by both, except for Kim 2 for which we 
use the more reliable measurement from Kim et ah (2015b) 
based on follow-up observations. We also use the method to 
predict Eri II’s PM, but caution that its large distance from 
the MW makes it unlikely that it is a satellite gravitationally 
bound to the MW. 

The resulting PM predictions are compiled in Table 4 
and illustrated in Figure 4. The table also provides the angle 
between the Galactocentric position of the object and the 
VPOS-3 plane, ^vpos -35 which is equal to the minimum 
possible inclination between the orbit of the object and the 
VPOS-3 plane. These angles tend to be small, again illus¬ 
trating that most of the new objects are closely aligned with 
the satellite plane. 

If future proper motion measurements reveal that more 
MW satellites follow the strong velocity alignment of the 
classical MW satellites with the VPOS, this would dra¬ 
matically increase the significance of their alignment. While 
the probability of a position drawn from an isotropic dis¬ 
tribution to be within an angle 0 from a given plane is 
Pvector = sin(^), the probability that a randomly oriented 
orbital plane is aligned to within 0 with a given plane is 
only Ppiane = 1 “ cos(^). For example, while 71 per cent of 
isotropically distributed positions are expected to lie within 
45° of a randomly oriented plane, only 29 per cent of isotrop¬ 
ically distributed orbital poles are expected to be within 45° 
of the plane’s normal vector. 


3.1 Interpreting PM predictions 

We caution against interpreting the PM predictions too 
tightly. The predicted PM ranges are derived from assum¬ 
ing the best possible alignment of the orbital pole with the 
normal to a satellite plane. Because the VPOS has a finite 
rms height, the orbits of most satellites can not be perfectly 
aligned. The orbital poles of the MW satellites therefore 
have a minimum intrinsic scatter around their average di¬ 
rection. In addition, all satellites which have orbital planes 
inclined from the VPOS will at some point during their or¬ 
bits be perfectly aligned with the VPOS. If observed at that 
position we would predict the satellites to orbit perfectly 
within the VPOS, even though they do not. 

The width of the VPOS can give us an idea of this in¬ 
trinsic scatter in orbital pole directions (to which scatter 
due to PM measurement uncertainties will be added). The 
VPOS-3 which is used for the PM prediction has a rms axis 
ratio of c/a = 0.2. This gives an approximate opening an¬ 
gle^ of 2 X arcsin(0.2) 23°. This estimate agrees with 
the measured spherical standard deviation of the eight best- 
aligned orbital poles derived from the observed PMs, which 
is about 27°. To meaningfully compare the predicted and 
observed PMs therefore requires that we measure the angle 
between the orbital poles (or the measured orbital pole and 
the VPOS normal vector assumed for the prediction). If the 


^ Strictly speaking the angle depends on (and decreases with) the 
radial distance from the MW, because the VPOS is a plane-like 
structure. 


inclination is less than the expected scatter, the measure¬ 
ment can be said to agree with the prediction. 


4 PREDICTED VELOCITY DISPERSIONS 

Here we attempt to predict the internal velocity dispersions 
of the newly discovered satellites from their reported photo¬ 
metric properties. Since it is unclear whether some of these 
objects are dwarf satellite galaxies or star clusters, we con¬ 
sider both possibilities. We consider both conventional grav¬ 
ity and MOND (Milgrom 1983). Since differences in the pho¬ 
tometric properties reported for the objects translate into 
differences in the predicted velocity dispersions, we make 
predictions using each of the reported sets of properties. 
Our results are compiled in Table 5 in the order in which 
they are listed in Table 1. 


4.1 Star Clusters 


If the newly discovered objects are star clusters, then they 
may be devoid of dark matter. In this case, their velocity 
dispersions follow directly from the virial relation and the 
observed stellar mass and half-light radii: 


a ~ 


3ri/2 J 


( 1 ) 


where ri /2 is the 3D half light radius (estimated as 4/3 the 
effective radius) and M* = To estimate the stellar 

mass, we assume a mass-to-light ratio of T* = 2 Mq/Lq. 
This is extremely uncertain (see §4.4.2). The velocity dis¬ 
persions predicted in this way are given in Table 5. These 
are low luminosity systems, so the predicted velocity disper¬ 
sions are small: < I kms“^ in all cases, and < 0.5 kms“^ 
in all but one case. However, the disruption of such dark 
matter free systems on their orbits around the Milky Way 
can increase their apparent mass-to-light ratios substantially 
(Kroupa 1997). This successfully predicted an object like the 
later discovered MW satellite Hercules (see the discussion in 
Kroupa et al. 2010). 


4.2 Scaling Relations 

If the newly discovered objects are dwarf satellite galaxies, 
then we expect them to reside within dark matter sub-halos. 
If this is the case, the kinematics are presumably dominated 
by dark matter as with the other known dwarfs. We then 
anticipate higher velocity dispersions. 

There is no universally agreed method to predict the 
velocity dispersions of individual dwarf satellite galaxies in 
ACDM. The correlation between luminosity and halo mass 
is exceedingly weak on these scales (Wolf et al. 2010), while 
the expected scatter in ACDM sub-halo properties is large 
(Tollerud et al. 2011). Consequently, the observed luminosity 
should have little predictive power, with the velocity disper¬ 
sions being essentially stochastic. 

There are empirical scaling relations that we can use 
to anticipate the velocity dispersion of the newly discovered 
objects. This is done with the scaling relations of McCaugh 
et al. (2007) and Walker et al. (2009), the results of which 
are tabulated in Table 5. Walker et al. (2010) showed that 
the distinct relations of McCaugh et al. (2007) and Walker 
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Table 5. Predicted velocity dispersions for each set of photometric parameters compiled in Table 1 


Name 

CNewton 

<^M07 

[km 

<7W09 

s-1] 

cmond 

c^MOND 

^Newton 

din 

[km^ s" 

C-ext 

kpc“^] 

Kim 1 

0.14 

1.4 

1.2 

0.18 

1.3 

110 

2280 

Ret II 

0.38 

3.2 

1.6 

0.68 

1.8 

120 

1160 


0.26 

3.0 

1.6 

0.46 

1.8 

90 

1210 

Tri II 

0.17 

3.1 

1.6 

0.32 

1.9 

60 

1060 

Tuc II 

0.24 

5.9 

2.1 

0.56 

2.4 

40 

660 


0.23 

7.6 

2.3 

0.61 

2.6 

30 

540 

Hor II 

0.21 

3.7 

1.7 

0.61 

2.9 

60 

430 

Hor I 

0.28 

4.2 

1.8 

0.86 

3.1 

70 

390 


0.38 

2.9 

1.6 

1.11 

2.9 

130 

430 

Phe II 

0.41 

3.1 

1.6 

1.30 

3.2 

140 

370 


0.30 

2.8 

1.5 

0.89 

3.0 

110 

430 

Eri III 

0.39 

1.8 

1.3 

1.27 

3.2 

230 

350 


0.25 

2.3 

1.4 

0.78 

3.1 

120 

390 

Kim 2 

0.24 

1.9 

1.3 

0.79 

3.3 

130 

350 


0.34 

1.9 

1.3 

0.89 

2.6 

190 

550 


0.34 

3.4 

1.7 

1.10 

3.2 

110 

360 

Gru I 

0.25 

4.5 

1.9 

0.88 

3.6 

60 

290 

Pic I 

0.36 

3.5 

1.7 

1.35 

3.7 

110 

270 


0.32 

3.0 

1.6 

1.15 

3.6 

110 

290 

Hyd II 

0.47 

4.4 

1.9 

1.78 

3.8 

110 

270 

Peg HI 

0.27 

4.7 

1.9 

1.30 

4.8 

70 

160 

Eri II^ 

1.00 

6.7 

2.2 

3.3 

3.3 

160 

100 


0.68 

7.0 

2.2 

2.8 

4.1 

100 

90 


crNewton^ Predicted velocity dispersion assuming M/L = 2 and that the objects are dark-matter-free. 

(PMor'- Predicted velocity dispersion assuming the empirical dark matter halo scaling relation of McGaugh et al. (2007). 
(Twog: Predicted velocity dispersion assuming the empirical dark matter halo scaling relation of Walker et al. (2009). 
<^mond: Predicted MONDian velocity dispersions, assuming stellar mass-to-light ratios of M/L = 2 for all satellite objects. 
<^MOND/(^Newton: Ratio of predicted MONDian to predicted Newtonian velocity dispersion. 

Uin: Internal MONDian acceleration of the satellite at its half-light radius, 
next: External acceleration acting on the satellite due to the potential of the MW. 


^ Eri II is the only object in the isolated MOND regime, i.e. for which next < «in- 


et al. (2009) were consistent with each other and the data 
available at the time. However, they are not identical, and 
give different results when extrapolated into the regime rep¬ 
resented by the new objects. We therefore tabulate distinct 
predictions for each, bearing in mind that these are the ex¬ 
trapolations of empirical scaling relations and are not pre¬ 
dictions derived from a specific theory like ACDM. 

McGaugh et al. (2007) fit the baryon subtracted rota¬ 
tion curves of spiral galaxies to obtain an expression for the 
rotation velocity due to the dark halo component. Though 
fit to spirals at larger radii, this relation did a good job^ 
of anticipating the common mass scale found for satellite 
galaxies by Strigari et al. (2008). Assuming a = VhlV^ and 
evaluating at the observed half light radius (in kpc), the 
relation of McGaugh et al. (2007) becomes 

log cr = 1.23 + 0.5 log re- (2) 


^ Evaluation of the scaling relation of McGaugh et al. (2007) at 
300 pc anticipates M(< 300 pc) = 1.8 x 10^ Mq. 


This anticipates velocity dispersions in the range of 1 - 
7 kms“^ (Table 5). 

Walker et al. (2009) fit the data for dwarf satellite galax¬ 
ies to obtain 


log cr 0.5 + 0.2 log re- (3) 

This should be more applicable to the newly discovered ob¬ 
jects, if they are indeed dwarf satellites, albeit very small 
ones. The dependence on size in this regime is rather weaker, 
anticipating velocity dispersions in the range 1 - 2.3 kms“^ 
(Table 5). Note that both scaling relations are known to be 
violated in some cases (Gollins et al. 2014). 

We may of course have a heterogeneous mix of objects: 
some might be dwarf satellite galaxies in sub-halos, while 
others might simply be star clusters. If so, the different an¬ 
ticipated velocity dispersions should help distinguish these 
two cases. However, considerable observational care will be 
required to do so given the small anticipated dispersions. 
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4.3 MOND 

We also tabulate the velocity dispersions predicted by 
MOND (Milgrom 1983). In this theory, the velocity dis¬ 
persion should follow from the observed properties of each 
object, irrespective of whether it is a star cluster or dwarf 
satellite galaxy. Either way, the physics is the same. 

To predict velocity dispersions with MOND, we follow 
the procedure outlined by McGaugh Sz Milgrom (2013a). 
We assume (as done there) that T* = 2 Mq/Lq. This 
approach has had considerable success in predicting, often 
a priori^ the velocity dispersions of the dwarf satellites of 
Andromeda and of isolated dwarfs in the LG (McGaugh Sz 
Milgrom 2013b; Pawlowski Sz McGaugh 2014b). A different 
method employed by Liighausen, Famaey Sz Kroupa (2014) 
has also been successful in predicting the velocity dispersions 
of the most luminous MW dSph satellites Fornax and Sculp¬ 
tor, but found the measured velocity dispersions of Sextans, 
Garina and Draco to be higher than predicted. 

The velocity dispersion predicted by MOND depends 
on whether the internal gravitational field ain of an object 
dominates (the isolated case), or if it is dominated by the 
external field effect (EFE) of the host galaxy aex- The case 
that applies depends on the relative strength of the inter¬ 
nal and external fields: a system is considered to be in the 
MOND regime and isolated if aex < CLin < ao, and is in 
the EFE regime if ain < CLex- There are thus two velocity 
dispersion estimators: that for the isolated case (eq. 2 of Mc¬ 
Gaugh & Milgrom 2013a), and that for the EFE case (eq. 3 
of McGaugh & Milgrom 2013a). 

The internal field depends only on the properties of each 
object, and is estimated at the half light radius as in Mc¬ 
Gaugh Sz Milgrom (2013a). The external field depends on 
the total baryonic mass of the Milky Way. For specificity, 
we adopt the empirical (“bumps Sz wiggles”) Milky Way 
model of McGaugh (2008), which provides an estimate of 
dex = yiiwlR cit the galactocentric distance R = dMW of 
each object. The uncertainty in the circular velocity of the 
Milky Way at these distances affects the predicted velocity 
dispersions at the ±0.3 kms“^ level. 

The EFE dominates in all but the most distant case: 
Eri II is the only object in Table 5 in the isolated regime. In 
some cases, the external field is only marginally dominant. 
In these cases, neither mass estimator is really adequate, 
and the velocity dispersion may be slightly under-predicted 
(Milgrom 1995). In all cases, the enhancement in the velocity 
dispersion predicted by MOND is rather modest. Typically 
it is only a factor of a few above the purely Newtonian (star 
cluster) case, and is often less than anticipated by the dark 
halo scaling relations. The many uncertainties of the quan¬ 
tities that go into the prediction are compounded by the 
potential for systematic errors. 


4.4 Challenges to Interpretation 

Gan we hope to observationally distinguish between the var¬ 
ious predictions? An accuracy of ~ 0.1 kms“^ is required to 
resolve the velocity dispersions anticipated for star clusters. 
This is certainly possible, if challenging. 

The larger concern is systematic uncertainties. For ex¬ 
ample, are these tiny systems in dynamical equilibrium? If 
not, we might misinterpret a high velocity dispersion of a dis¬ 


solving star cluster as the equilibrium dispersion of a dwarf 
satellite residing in a dark matter sub-halo. 

Many of the uncertainties involved in measuring and in¬ 
terpreting velocity dispersions have been discussed by Mc¬ 
Gaugh Sz Wolf (2010). We emphasize here just two possi¬ 
ble systematics that we fear will make it extremely difficult 
to distinguish between the various possibilities. One, binary 
stars may inflate the observed velocity dispersions. Two, the 
conversion from light to stellar mass is rather fraught for 
systems containing so few stars. 

4 . 4.1 Binary Stars 

Some of the individual stars for which velocities are ob¬ 
tained are presumably members of unresolved binaries. Bi¬ 
nary stars are themselves in orbit around one another, possi¬ 
bly at speeds comparable to the velocity dispersion expected 
for the system as a whole. This can inflate the velocity dis¬ 
persion measured for the system. Binaries can easily con¬ 
tribute enough to the measured velocity dispersion to change 
the interpretation from one extreme to the other. One can 
correct for this effect, but it requires the patience of many 
repeat observations (Simon et al. 2011). 

4.4- 2 Stellar Mass 

Another problem when considering ultrafaint dwarfs is the 
uncertainty in the stellar mass-to-light ratio. Ultrafaint 
dwarfs are so small that their entire luminosity can be less 
than that of a single high mass star. This violates an es¬ 
sential assumption in the estimation of stellar mass-to-light 
ratios with stellar population models: that there are enough 
stars to statistically sample all phases of stellar evolution. 
Indeed, in systems composed of only a few hundred stars, 
the evolution of a single star up the giant branch will sub¬ 
stantially change the luminosity of the entire system without 
changing its mass. For this reason, all predictions made with 
an assumed mass-to-light ratio are subject to large uncer¬ 
tainty. This strongly affects the predictions for both purely 
Newtonian star clusters and MOND. 

4.4- S Ret II, Hor I and Hyd II 

Ret II is the first of the objects in Table 5 to have a 
measured velocity dispersion a and heliocentric velocity 
^hei- Walker et al. (2015) measure a — d.b^Q'g kms“^ 
(^hei = 64.8lJ‘Jkms“^), Simon et al. (2015) measure 
a — 3.3 ± 0.7 kms“^ (uhei = 62.8 ± 0.5kms“^) and Ko- 
posov et al. (2015b) measure a = 3.2lJ‘5 kms“^ (uhei = 
64.7lj g kms“^). These measurements are nicely consistent, 
and are clearly too large for a star cluster devoid of dark mat¬ 
ter (0.3 kms“^) or for MOND (0.5 kms“^). This of course 
presumes that the system is in dynamical equilibrium, that 
our guess for the stellar mass-to-light ratio is not far off, and 
that binary stars contribute ^ 3 kms“^ in quadrature to 
the observed velocity dispersion. 

We have used the opportunity provided by the indepen¬ 
dent measurements of Ret IBs velocity dispersion to check 
whether the studies are in mutual agreement or whether sys¬ 
tematic errors might be present. This is motivated by the 
comparison for the Garina dSph by Godwin Sz Lynden-Bell 
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Figure 5. Velocity measurements for those Ret II stars in com¬ 
mon between the sets Walker et al. (2015), plotted on the vertical 
axis, and Koposov et al. (2015b), plotted on the horizontal axis. 
The data points are plotted with their respective error bars, also 
shown is the best fitting line and the upper and lower 2-sigma 
confidence bands. We find that the measurements agree well with 
each other, giving confidence in the reported heliocentric velocity 
and velocity dispersion for this object. 


(1987). We can make a similar comparison for the studies of 
Walker et al. (2015) and Koposov et al. (2015b), who have 
13 Ret II stars in common. The respective velocities mea¬ 
sured for these stars are plotted against each other in Fig. 
5. We find that the best-fitting line has a slope of 0.79 zb 0.21 
and a y-axis intercept of 13.6zbl3.3 (not counting the outlier 
with the largest error bars which also disagrees with the sys¬ 
tematic velocity of Ret II by almost 15 kms“^ in one of the 
two studies). Hence, there is indeed a positive correlation 
between the two velocity sets, which furthermore is consis¬ 
tent with a slope of one. This supports the interpretation 
that not random errors but indeed the internal velocity dis¬ 
persion of the object is a major contributor to the spread in 
velocities. 

Taken at face value, the observed velocity dispersion is 
consistent with the empirical scaling relation of McGaugh 
et al. (2007), which anticipates cr = 3 - 4 kms“^, depend¬ 
ing on whose photometry is employed. This implies a dark 
matter halo consistent with the near-universal halo found 
by McGaugh et al. (2007) and Walker et al. (2010). Unfor¬ 
tunately, this empirical dark matter halo is not consistent 
with AGDM (McGaugh et al. 2007), though presumably it 
can be accommodated by invoking feedback or some other 
mechanism. 

To make matters worse, the universal halo that is suc¬ 
cessful in the case of Ret II does not work in the cases of the 
satellites of M31 And XIX, XXI, and XXV (Gollins et al. 
2014). These objects are faint, but have much larger effec¬ 
tive radii than the objects under consideration here. Appli¬ 
cation of equation 2 of McGaugh et al. (2007) anticipates 
cr > 13 kms“^ for these dwarfs of M31, while they are ob¬ 


served to have cr < 5 kms“^ (Gollins et al. 2014). These 
objects should be strongly affected by the EFE in MOND, 
which was unique in accurately predicting their velocity dis¬ 
persions in advance (McGaugh Sz Milgrom 2013a,b). We 
therefore urge caution in interpreting the velocity disper¬ 
sions of these objects, especially in light of the systematic 
uncertainties discussed above. 

Koposov et al. (2015b) also infer a velocity dispersion 
for Horl from five stars of cr = 4.91q'9 kms“^ (r’hei = 
112.8l2'6 knis“^). This dispersion value exceeds all predic¬ 
tions, but it again comes closest to the prediction using the 
McGaugh et al. (2007) scaling relation. However, for the 
reasons discussed above one should be extremely cautious 
in interpreting these velocity dispersion measurements. 

After this manuscript was submitted, Kirby, Simon & 
Gohen (2015) announced the first spectroscopic measure¬ 
ment of stars in Hydll. They did not resolve its velocity 
dispersion, but report an upper limit of cr < 4.5kms“^ 
(95per cent confidence), which is consistent with all pre¬ 
dictions. Kirby, Simon & Gohen (2015) found a heliocentric 
velocity for Hydll of Vhei = 303.1 zb 1.4kms“^, which they 
report to be similar to the Leading arm of the Magellanic 
Stream. Such a similarity to the velocity of the Magellanic 
Stream is a general trend for the objects found to lie within 
the dwarf galaxy planes in the LG (Pawlowski, Kroupa Sz 
Jerjen 2013), which is also followed by Retll and Horl (see 
Eig. 6). Unfortunately, the line-of-sight velocity of satellite 
objects does not provide decisive information on whether it 
orbits within the VPOS, because this velocity is mostly ori¬ 
ented along the radial component of its position vector from 
the Galactic center. 


5 DISCUSSION AND CONCLUSIONS 

We have compiled a list of 14 recently discovered stellar sys¬ 
tems in the vicinity of the MW, many of which are probably 
MW satellite galaxies. We find that most of these objects 
align well with the VPOS, which consists of both satellite 
galaxies and star clusters (Pawlowski, Pflamm-Altenburg Sz 
Kroupa 2012). The updated VPOS fit parameters, compiled 
in Table 2, do not deviate substantially from the previous 
ones: the rms height is almost unchanged, the orientation 
preserved to within 9° and the offset from the MW center is 
reduced. Assuming that this alignment indicates the objects 
to be part of a common dynamical structure, as is indicated 
by the aligned orbital poles of the 11 classical MW satellites 
(Pawlowski Sz Kroupa 2013), we predict the proper motions 
of the new satellite objects (see Sect. 3). 

We apply Newtonian and MONDian dynamics and dif¬ 
ferent dark matter halo scaling relations to predict the veloc¬ 
ity dispersions of the objects from their photometric proper¬ 
ties (see Sect. 4). These three distinct assumptions result in 
predictions with only modest differences. Eor most objects 
Newtonian dynamics predict velocity dispersion between 0.2 
and 0.4kms“^, the dark matter scaling relations predict ve¬ 
locity dispersions between 1 and 4kms“^, and the MOND 
predictions lie in between these. This small range of very 
low velocity dispersions makes it extremely difficult to dis¬ 
criminate between the three cases observationally, which will 
require very precise measurement and good control of sys¬ 
tematic effects such as unresolved binary stars. 
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Figure 6. Comparison of the Magellanic Stream (from Nidever et al. 2010) and the LG dwarf galaxies as in figure 17 of Pawlowski, 
Kroupa & Jerjen (2013), but updated to include the newly discovered objects and additional information. The upper panel plots positions 
on the sky in the Magellanic Stream Coordinate (Nidever, Majewski &: Burton 2008). The solid blue line indicates the intersection of the 
best-fit VPOS plane with a sphere of 250 kpc radius. The arrows indicate the current most-likely velocity vectors of the MW satellites 
LMC, SMC, Fornax, Carina and Sculptor, as compiled in table 2 of Pawlowski & Kroupa (2013). The velocity vectors are represented 
in position space by converting 5kms“^ to length of 1 kpc and then projected into Magellanic Stream coordinates. As in Fig. 1 the 
velocity vectors reveal that the satellites move mostly along the plane. The blue plus signs indicate the positions of the MW satellites 
(dark blue for the two star clusters Kim 1 and 2), yellow dots the positions of the non-satellite dwarf galaxies which are part of a highly 
flattened plane approximately connecting the MW and M31, and crosses are satellite galaxies of M31 (red for those in the M31 satellite 
galaxy plane). The black ellipse indicates the position and orientation of M31. The lower panel gives the line-of-sight velocities of the 
Magellanic Stream and those galaxies for which kinematics are available. See Pawlowski, Kroupa &: Jerjen (2013) for further details. 


Most of the objects are in the southern hemisphere of 
the MW, and the majority of them have been discovered in 
the DES. The survey footprint lies close to the Magellanic 
clouds which orbit within the VPOS, such that an align¬ 
ment with the VPOS might not be unexpected. However, 
even though the area covered by the DES so far falls close 
to the VPOS, one would expect the ten objects discovered in 
the data to have about 50 per cent larger mean and median 
offsets from the VPOS if they were drawn from an isotropic 
distribution confined to the survey footprint. Both mean 
and median offsets at least as small as observed are rare 
among such randomized realisations (9 and 4 per cent, re¬ 
spectively). Eurthermore, several other objects were discov¬ 
ered elsewhere, away from known concentrations of satellite 
galaxies, but nevertheless aligned with the VPOS. The fact 
that the PanSTARRs survey, despite its Stt sky coverage, 
has so far not resulted in the discovery of a large number 
satellite galaxies outside of the VPOS provides further hints 


that even the fainter MW satellites align with the satellite 
structure which was first discussed almost 40 years ago by 
Kunkel & Demers (1976) and Lynden-Bell (1976). 

Among the M31 satellite galaxies a similar, and appar¬ 
ently also co-rotating, plane consisting of about half of the 
satellite population was found by Ibata et al. (2013). It is 
aligned with the M31’s prominent stellar streams (Hammer 
et al. 2013), which is reminiscent of the preferential align¬ 
ment of streams in the MW halo with the VPOS, most 
prominently the Magellanic Stream. Co-orbiting planes and 
similar satellite alignments might even be common through¬ 
out the Universe (Ibata et al. 2014; Pawlowski & Kroupa 
2014; Tully et al. 2015). 

The new satellites could prove to be important for the 
wider picture of dwarf galaxies and the overall dynamics in 
the Local Group (LG). The non-satellite LG dwarf galaxies 
are confined to two highly symmetric and extremely narrow 
planes (Pawlowski, Kroupa & Jerjen 2013). The dominant 
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of these two planes appears to connect M31 and its satellite 
galaxy plane with the VPOS around the MW, and agrees in 
projected position and line-of-sight velocity with the Magel¬ 
lanic Stream (see section 7.4 in Pawlowski, Kroupa Sz Jerjen 
2013 and Fig. 6). The new discoveries are particularly inter¬ 
esting because many lie in the vicinity of the Magellanic 
Clouds but into the direction of M31, close to the Magel¬ 
lanic Stream. This region was identified as the “direction of 
decision” by Pawlowski, Kroupa & Jerjen (2013), because 
it is where the MW and M31 satellite planes intersect with 
the dominant plane of non-satellite dwarf galaxies in the LG. 
Knowledge of the phase-space distribution of objects in this 
region should help to determine if and how these structures 
are connected. 
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